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ABSTRACT 

To provide a quantitative cosmological context to ongoing observational work on the 
formation histories and location of compact massive galaxies, we locate and study a 
sample of exceptionally compact systems in the BOLSHOI simulation, using the dark 
matter structural parameters from a real, compact massive galaxy (NGC 1277) as 
a basis for our working criteria. We find that over 80% of objects in this nominal 
compact category are substructures of more massive groups or clusters, and that the 
probability of a given massive substructure being this compact increases significantly 
with the mass of the host structure; rising to ^ 30% for the most massive clusters in 
the simulation. Tracking the main progenitors of this subsample back to z = 2, we 
find them all to be distinct structures with scale radii and densities representative of 
the population as a whole at this epoch. What does characterise their histories, in 
addition to mostly becoming substructures, is that they have almost all experienced 
below-average mass accretion since z = 2; a third of them barely retaining, or even 
losing mass during the intervening 10 Gyr. 
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1 INTRODUCTION 

Massive (M* .li, 10^^ M©) compact (Ro%1.5kpc) galaxies 
are extremely rare in the present-day Universe, several 
groups estimating their number density to be currently 
n < 10“®Mpc“^ (Trujillo et al. 2009; Taylor et al. 2010). 
As well as for their scarcity, these objects are of great inter¬ 
est because their structural properties (Trujillo et al. 2012) 
are comparable to a large fraction of the massive galaxy pop¬ 
ulation 10 Gyr ago at a redshift of 2 ~ 2 (e.g. Trujillo et al. 
2007; Buitrago et al. 2008; Carrasco et al. 2010; Chevance 
et al. 2012; Szomoru et al. 2012), when compact massive 
galaxies were significantly more abundant ( > 10“'*Mpc“^; 
Taylor et al. 2010). 

As cosmic time has evolved, these galaxies have progres¬ 
sively disappeared (Stockton et al. 2010; Valentinuzzi et al. 
2010; Cassata et al. 2011; Poggianti et al. 2013; Damjanov 
et al. 2014; Hsu et al. 2014). The decline in the number of 
these galaxies has been linked with merging (Lopez-Sanjuan 
et al. 2012; Newman et al. 2012; Oser et al. 2012; Bedorf 
& Portegies Zwart 2013; Hilz et al. 2013; Quilis & Trujillo 
2013; Shankar et al. 2013) and this scenario can explain a 
large number of observations (e.g. Trujillo et al. 2011, and 
references therein). 

The recent discovery of a nearby compact massive ‘relic’ 
© 2015 RAS 


galaxy, NGC 1277, with similar structural properties to typ¬ 
ical galaxies of the same mass at z~2 (Trujillo et al. 2014), 
has opened the possibility of exploring the properties of this 
galaxy population in unprecedented detail, including in par¬ 
ticular some hints about the dark matter halo properties 
(van den Bosch et al. 2012). Having access to these char¬ 
acteristics is one key to understanding this type of object 
in a cosmological context. In fact, we lack knowledge about 
the connection of these compact massive galaxies with the 
cosmic structures which they inhabit. In particular, which 
kind of dark matter halo they occupy and how these halos 
evolve with time. 

In this paper we aim to fill this gap, and address the 
following questions: How do these particular dark matter 
structure properties compare with the rest of the popula¬ 
tion? Are they substructures of larger dark matter halos or 
are they isolated objects? How have the structural proper¬ 
ties of these dark matter halos evolved with time? 

We confront these questions from a cosmological per¬ 
spective by taking the dark matter structure properties that 
would be consistent with our knowledge of this real system 
(§2), and then looking for structures in a numerical simu¬ 
lation that share these properties (§3.1). These candidate 
structures thus identified, we can then begin asking where 
they came from and what has dictated their evolution (§3.2). 
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2 ASSOCIATING COMPACT GALAXIES AND 

COMPACT STRUCTURES 

We begin by characterising the properties of one particularly 
well-known compact massive galaxy, NGC 1277, described 
in Trujillo et al. (2014) to have a stellar mass of M* = 1.2 ± 
0.4 X 10^^ Mq, half of which is contained within a circularised 
radius of just = 1.2 ± 0.1 kpc. 

To make the association between real systems such as 
this, and the structures in a cold dark matter (CDM) sim¬ 
ulation, we can draw on both empirical, semi-empirical and 
theoretical arguments. 

Ideally, we might rely entirely on the empirical mea¬ 
surement of the real system’s mass profile, but unfortu¬ 
nately the profile is only marginally constrained by direct 
measurement. However, we can still note the best fitting 
profile parameters, found from a full analysis of the stel¬ 
lar kinematics by van den Bosch et al. (2012). These are: 
pB ~ 0.027 Mqpc”^ and Cg = 26 kpc, where these apply to 
the ‘NFW profile’, the standard fitting function of Navarro 
et al. (1995, 1996, 1997) for density as a function of radius 
(Al). 

To use these figures directly to seek analogue structures 
in CDM-only simulations is of course to neglect the effects 
of the formation of the galaxy itself on the profile. But given 
the already stated uncertainty in ps and Vs, a detailed at¬ 
tempt to analytically reverse-model the effects of collapse 
and ejection of the normal matter component would be in¬ 
appropriate. Appealing also to numerical work on this issue, 
a recent study by Schaller et al. (2014) compared the density 
prohles of structures in the hydrodynamical EAGLE simula¬ 
tions (Schaye et al. 2015) with the equivalent structures from 
a ‘dark matter only’ version. Though both Vs and ps did dif¬ 
fer between the two cases, with standard deviations at the 
relevant mass range of about a factor of two, the expectation 
change was always close to zero. There was no clear tendency 
for the absence, or presence, of normal matter in the numer¬ 
ical simulations to shift these parameters systematically one 
way or other. 

To add to these empirical measurements of profile pa¬ 
rameters, we can also use the stellar mass of NGC 1277 to 
indicate a likely structure mass range, by appealing to the 
abundance-matching hypothesis (e.g. Behroozi et al. 2010; 
Moster et al. 2010). In particular, the expected halo mass 
for a given stellar mass^, calculated in this way by Behroozi 
et al. (2010), would associate our stellar mass range here, 
M* « 0.8 — 1.6 X IO^^Mq, with halo masses of: 

Mv « 2 X 10^^ - IO^^Mq , (1) 

where this approximate range incorporates both the uncer¬ 
tainty in M* and the scatter in the mapping to Mv. 

As well as using the stellar mass as a guide, we can also 
draw on standard theoretical arguments which would ex¬ 
pect the radial extent of the final galaxy to broadly reflect 
the quantity and distribution of angular momentum in the 
structure in which it formed (e.g. Fall & Efstathiou 1980; 
Fall 1983; Mo et al. 1998). This idea has recently received 
a boost of empirical support from Kravtsov (2013), who 
showed that galactic radii and host structure radii linked by 
the abundance-matching ansatz are directly proportional: 

^ as opposed to the expected stellar mass for a given halo mass. 


(Ro) ~ 0.015 (R 200 ) with a scatter of « 0.2 dex (R 200 being 
the radius that encloses 200 times the critical density). 

This correlation would pair a galaxy of Re ~ 1.2 kpc 
with a structure of R 200 ~ 80 kpc, which at 2 ; = 0 corre¬ 
sponds to M 200 ~ 200Ro^2oo/ 2G « 6 X lO^^M©, or a virial 
mass of approximately AR ~ lO^^M©. The discrepancy be¬ 
tween this mass, and the range in the preceding paragraph 
deduced from the stellar mass, is a reflection of the excep¬ 
tional properties of NGG 1277. 

To reconcile this, we need to go beyond the theory’s 
application to mean or representative values. The standard, 
more detailed version is that the net specific angular momen¬ 
tum of material is conserved in the formation of the galaxy. 
This would mean that a more compact halo, with more of its 
mass at smaller radii, leads to a galaxy formed from mate¬ 
rial with less angular momentum - and thus more compact 
- than others of the same mass. 

For the purposes of our exercise here, it would therefore 
seem sufficient to consider a category of compact structures 
whose mass exceeds the lower limit of the range indicated in 
(1), have highly centrally concentrated mass distributions, 
but with limits chosen such that our real system would be 
comfortably included in the sample (based at least on the 
best fitting parameters which are available). 

We thus define our sample of interest in terms of char¬ 
acteristic^ or inner density, pa, and total virial mass, AR: 

Pa > 0.02 Mqpc”^ and A/v > 2 x 10^^ M© , (2) 

This total inner density value corresponds to a dark matter 
density of 0.02(1 — flb/flM) ~ 0.016, so the best-fitting dark 
matter inner density for NGG 1277 would fall within this 
category by about a factor of two. The sample also contains 
the best fitting scalelength by a similar margin, as can be 
seen from the upper right panel of Fig. 1. In terms of con¬ 
centration parameter, c = r^/rs, the criteria in (2) amounts 
to a selection of about c R 20. 

So, based on all these combined considerations, we have 
a well-motivated and interesting category of compact, mas¬ 
sive structures. The next step is to locate such structures 
in a large numerical simulation and investigate their where¬ 
abouts and histories. This now follows in §3. 


3 ANALOGUE COMPACT STRUCTURES IN 
SIMULATIONS 

The following results and figures are taken from the publicly 
available catalogues of structures found within the BOLSHOI 
simulation (Klypin et al. 2011; Behroozi et al. 2013b, and ref¬ 
erences therein) using the halo finding algorithm presented 
in Behroozi et al. (2013a). The simulation is a cubic box of 
250h“^ comoving Mpc containing 2048® particles of mass 
just under 2 x 10® Mq . 

Cosmological parameters were chosen to be consistent 
with the WMAP five year (Dunkley et al. 2009; Hinshaw 

^ ‘Characteristic density’ is the term used by Navarro et al. 
(1997). The mean enclosed density for their profile (Al) falls to 
this value, pa, at a radius r = 0.72 rg. As ^ 5 — 10rg, this 
density pa is representative of a fairly centralised region, typically 
just a few percent of the total volume. Hence also ‘inner density’ 
(Bullock et al. 2001; Wechsler et al. 2002, and here). 
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Figure 1. The position of compact structures relative to the overall distributions in the BOLSHOI simulation, and to the position of 
their host structures, where relevant. The top panels show the the Mv — Tg plane (left) and the pg — rg plane (right) using the same 
key. Our category of ‘compacts’ are highlighted as dots, their hosts in circles, and the distribution of the population is indicated with 
background shading. The mass and density defining our nominal selection criteria (2) and other relevant loci, as labelled, are highlighted 
with dot-dashed lines. The star in the upper right panel indicates the best-fitting profile parameters for NGC 1277. The bottom panels 
show the projected distributions on to the respective x-axes. 


et al. 2009; Komatsu et al. 2009) and seven year data 
(Jarosik et al. 2011), notably: hO = 0.7, = 0.27, ilh = 

0.0469, as = 0.82 and Us = 0.95. 

The virial masses, Mv, and radii, r^, for structures in 
the simulation are defined by Klypin et al. (2011) relative 
to the mean universal matter density at that redshift: 

GMv = V2Av(z)//o"nM(l + zfrl , (3) 

where Av( 2 ;) is the virial overdensity. For the adopted cos¬ 
mology, Av(0) = 360, which is equivalent to nMAv(O) = 97 
times the critical density. The nominal resolution limit of 
Mv « lO'" Mq (Behroozi et al. 2013a), is comfortably ex¬ 
ceeded by all the structures that we are concerned with in 
this study. 


3.1 Where are they now? 

The final distribution of masses, radii and inner densities of 
structures in the simulation is shown in Fig. 1, highlighting 
in particular the location of the 346 structures in the vol¬ 
ume (n ~ 7 X 10“®Mpc“^) that meet our nominal compact 
criteria defined by (2). This subsample constitutes the most 
compact 1% of all structures in the same mass range, as can 


be appreciated by their location at the very bottom edge of 
the Mv — Ta plane, shown in the upper left panel. 

The upper right panel shows the Ps — ra plane. The best- 
htting parameters for NGC 1277, indicated with a star, show 
that these values would place it inside our compact subsam¬ 
ple (though the limits were motivated by a combination of 
factors; see §2). 

The lower left panel of Fig. 1 shows the total mass dis¬ 
tribution of the population and the subsample, divided by 
distinct systems and substructures. This provides a key re¬ 
sult that over 80% of such compact massive systems are 
identihed to be substructures of a larger collapsed region (as 
compared to 8% of all structures in the same mass range). 
Furthermore, of the ~ 20% of these exceptional systems that 
are in fact distinct, we note that over half are found to have 
passed through larger structures in the past. So, in total, 
94% of the objects in this subsample are or have been influ¬ 
enced by a larger collapsed region. 

This is in agreement with the overall results from this 
simulation, which showed a trend for substructures to be 
more compact on average than distinct systems of the same 
mass (Klypin et al. 2011, eqns. 10-11), and with earlier work 
(e.g. Ghigna et al. 1998; Bullock et al. 2001). It is also in 
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Figure 2. The distribution of the compact structures in our nom¬ 
inal definition (highlighted in Fig. 3) as a function of the mass of 
their host structure at 2 : = 0. Specifically, for each host structure 
we find the number of compact substructures within it and divide 
by the number of substructures that lie in the same mass range. 
Points show the mean of this fraction for all structures in each 
host mass bin that contain one or more substructures in the mass 
range. Outer, dashed error bars give the standard deviation in 
this fraction, and inner error bars indicate the error on the mean 
(i.e. a!y/n, where n is the number of host systems in the mass 
interval which contain substructures from the given mass range). 
Systems containing no substructures in the mass range are not 
included in the statistics. 

agreement with the actual location of our case study, NGC 
1277, in the massive Perseus Cluster. 

The hosts of these compact satellites^ are also high¬ 
lighted in Fig. 1, and are shown to have very flat distribu¬ 
tion in mass from 10^^ — M©. So given that you have a 

compact system it is similarly likely to turn out to be in any 
sufficiently massive host. 

The clause here is emphasised because the converse 
statement, is not true. Taking a potential host structure at 
random, it is overwhelmingly more likely to find one of these 
compact substructures if the host is more massive. This is 
clear when the number of hosts-of-compacts (hatched area in 
Fig. 1) is considered as a fraction of the number of systems 
in general (solid fine above it), which is rapidly declining. 
However this is to some extent a trivial result in the sense 
that there are just more objects inside the most massive clus¬ 
ters. Whatever it is you are looking for, it therefore seems 
a good bet that you are more likely to find it there, and 
compacts are apparently no exception. 

But with this large and detailed catalogue at our dis¬ 
posal, it should be possible to get a little beyond these two 
basic statistical effects and see if - after these have been 
accounted for - there remains any other evidence of a pre¬ 
ferred environment for this category of system (i.e. one that 
is not also the case for systems in general). 

For those that are still substructures, we seek some nor¬ 
malisation of probability that gives us a fair sense of their 
likely whereabouts. One such normalisation is to look at 
their abundance as compared with other substructures that 
fie in the same mass range as they do. This comparison is 
shown in Fig. 2. 

The result of this calculation is this: if one finds a sub- 

® About 10% of the hosts contain two of the compact systems 
(none contain three), so there are slightly fewer hosts than there 
are compact substructures. 


Figure 3. The distribution of systems in the mass—density plane 
(left panel), showing also the projection on to the density axis 
(right panel). Solid and dashed lines show the mean and standard 
deviation of the inner density of all structures in each mass bin. 
The dotted line shows the mean for substructures only. 

structure in the given mass range (~ 2 x 10^^ — 10^® M©), 
then the probability of it being from this highly compact cat¬ 
egory - the most compact 1% of all structures in that mass 
range - increases from a few to 30% as a rising function of 
the mass of the host structure. 

There is considerable scatter in the results, as indicated 
by the dashed error bars in Fig. 2. This is perhaps to be ex¬ 
pected given the low number statistics involved in the frac¬ 
tion itself; the number of such massive substructures con¬ 
tained in any one group or cluster is at most eight. Many 
hosts in the sample have one or more massive substructures 
but no ‘compacts’ (i.e. 0%) and some have just one massive 
substructure which is compact (i.e. ‘100%’). 

However, as there are over 2,600 such hosts within the 
simulated sample, the trend seen in Fig. 2 is statistically 
significant. This is reinforced by the solid error bars, which 
show the much smaller error on the mean, af ^Jn. Whilst 
a randomly selected structure from this mass range is only 
~ 1% likely to be so compact, a substructure from the same 
mass range in a 10^"* Mq cluster has about a 10% chance, 
rising to 20 — 40% in a lO'^® Mq cluster. 

These conclusions should be followed by reiterating that 
these figures are specific to our choice of category, that the 
percentages given do not include host structures which have 
no massive substructures, and that 20% of this compact cat¬ 
egory are in fact not substructures at all. But, that said. 
Figs. 1 and 2 reveal a great deal about these systems’ where¬ 
abouts in the single, z = 0, halo catalogue. 

Before moving on to higher redshifts, there are a few 
small points of interest concerning features of the general 
distribution of structures revealed in Fig. 1. Firstly, the 
mean scale radius as a function of mass, in the upper left 
panels, is not far from tracking a constant mean density so¬ 
lution, but is a little steeper reflecting the fact that, at low 
redshifts, less massive halos tend to be more concentrated 
(Klypin et al. 2011, fig. 5). A specific consequence of this 
shows up in the distribution of inner densities, to the right, 
showing that all the highest inner densities are found in 
small, lower mass systems. A glance at the lower right panel 
confirms also that all these very high inner density systems 
are substructures. 

Having mentioned these features, it is worth taking a 
quick look at the mass-density plane itself, shown here for 
reference in Fig. 3. This firstly reconfirms the simple point 
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Figure 4. The distributions of the structure population a.t z = 2, using the same axes and key as Fig. 1. The highlighted subsample 
here are the most massive progenitors of the systems that were highlighted at 2 = 0. 


that all the top three orders of magnitude in inner density 
(e.g. Xi 0.1 Mq/pc®) occur in low-mass substructures (i.e. 

A 10^^ M©), or former substructures. It also shows the more 
general feature of decreasing inner density with total struc¬ 
ture mass. The mean varies by a factor of ten for systems in 
general, with a signihcantly more pronounced variation of a 
factor of a hundred seen in substructures. 

Fig. 3 also includes a projection of distribution on to 
the density axis, as did Fig. 1, but in this case using a lin¬ 
ear scale. This emphasises that, notwithstanding these im¬ 
portant variations, the general distribution is quite sharply 
peaked and most inner densities lie within an order of mag¬ 
nitude of this overall mode. 

3.2 Where did they come from? 

Having located and characterised our subsample of compact 
massive structures at 2 = 0, it is immediately interesting to 
ask what the basic features of their histories are. This should 
help explain what has led to them being such rare objects, 
or at least assess any expectations we may have had already 
from basic structure theory. 

As a first step, we can begin by taking the same axes 
that were used to characterise the final properties of the 
structure population in Fig. 1, but show the state of affairs 
10 Gyr previously, at 2 = 2. This is shown in Fig. 4. 

Regarding first the two general distributions in these 


planes, these do not change a great deal from 2 = 2 —> 0. 
There is, of course, continued growth of high mass structure 
but these new systems fall on the existing trends with scale 
radius and density. So the overall distributions stay the same 
in that sense. 

The extremes of the density distribution become in¬ 
creasingly populated with time. The low-density end should 
not be over interpreted as the NFW fits are no longer re¬ 
liable as Ta —>■ Cv (see appendix A). The high-density end, 
though, is interesting. That this tail extends as time goes on 
is perhaps not as expected, and not all these extreme objects 
at 2 = 0 are simply older structures that have been stuck in 
a cluster and not evolving; many become more dense from 
2 = 2 —>• 0 (see also Fig. 6 and Fig. Bl). 

This is a good point to revisit our expectations concern¬ 
ing the relationship between the ‘age’ of structures and their 
density. The familiar statement of this relationship is that 
structures’ “characteristic densities are just proportional to 
the cosmic density at the time they ‘formed’ ” (Navarro et al. 
1997), where formation is taken to be the point at which 
“half of the final mass is in collapsed progenitors more mas¬ 
sive than 10% of the final mass”. Here, in Fig. 5, we take 
the opportunity to compare this statement with the much 
more recent simulation results we now have to hand. 

Despite using the slightly later time that is recorded 
in the BOLSHOI catalogue (the point at which half the peak 
mass lies in the most massive progenitor), the correlation 
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found in Fig. 5 does follow this well-known rule, at least for 
typical formation redshifts. This can be seen by comparing 
the two solid lines in the bottom panel; one representing 
the mean half-mass redshift of structures in the simulation, 
and the other showing the locus ps oc (1 -I- z)^, found by 
Navarro et al. (1997), normalised to go through the overall 
mean values here. The two lines are overlapping from 2 « 
2 —>■ 0.5, during which period the majority (« 60%) of all 
these Mv > 10^^ Mq structures attain half their peak mass. 

The mean half-mass redshift for all the structures in our 
final z = 0 sample turns out to be 2 ; ~ 1.3. The mean for 
the compact subsample is somewhat earlier than that, but 
this is not necessarily what is making them exceptional. The 
distribution of formation times is broadly representative of 
the population, so we must assume that it is the later half 
of their evolution which makes them exceptional by 2 = 0. 

With this in mind we now turn back to the compact 
candidates’ locations in Fig. 4, recalling that, at 2 = 0, these 
all end up in a very sparsely populated region at the edge 
of the distribution in Fig. 1. At 2 = 2, however, it seems 
that their most massive progenitors were evenly distributed 
around the mean trends in both planes. In fact, none of the 
progenitors fall within the region on the My — fs plane which 
their descendants will eventually occupy, and very few fall 
in this region in the pa — Ta plane. 

The only thing which these main progenitors of this 
subsample seem to have in common at all, at this earlier 
time, is their mass. Though the spread is considerable (over 
an order of magnitude) it is not actually that much greater 
than the range of their descendants at 2 = 0. And, crucially, 
the centres of both distributions are almost the same, mean¬ 
ing that the thing these compacts seem to have in common 
is very little mass growth since 2 = 2. 

To explore and quantify this, we can look specifically 
at the distribution of structures in terms of this property - 
mass growth - and see how this correlates with changes in 
inner density. This is shown in Fig. 6. 

Concerning first the properties of the population in gen¬ 
eral, this figure shows a clear anti-correlation between mass 
change and inner density change. As inner density for a given 
mass increases monotonically with concentration (see Ap¬ 
pendix B), Fig. 6 can be interpreted also in this sense; lesser 
mass growth tending to lead to a more concentrated final 
structure, and mass loss even more so. 

Fig 6 also shows that the modal'* density change is close 
to zero. This runs alongside the prominent idea of inside-out 
structure growth (Diemand et al. 2007; Cuesta et al. 2008), 
which can be further corroborated by looking also at the 
changing density inside a fixed physical radius, shown for 
reference in Appendix B. 

But, in particular when mentioning that notion, the 
scatter in all directions in this figure cannot be emphasised 
enough. All four quadrants of the main panel are occupied. 
It is fair to say that the true fraction of the population lying 
away from the origin is exaggerated somewhat by the log¬ 
arithmic scales, but we set this visualisation issue aside by 
including the simple statistic in the right panel that a clear 



10’’ 

V. 

p 

10“ 

y 

B 

1000 

0 

100 

6 

z 

ro 

1 

y 

n. 

10 

1 

© 

s 

10 

o' 

II 

1 


0.1 

c/3 

G 

y 

t3 

0.01 

u 

y 

G 

■B 

10“” 

G 

'iZ 

10““ 


Redshift at half peak mass, 2 , /, 

6 5 4 3 2 1 0 



_ (/?g) ±a-, all structures 


0.2 


0.4 0.6 0.8 1 


Scale factor at half peak mass, 


Figure 5. The correlation between the final (2 = 0) inner density 
of structures, and the redshift or scale factor at which half their 
peak mass was in place. The key is as labelled, following the theme 
set by Figs. 1-6, with the addition of a faint solid line showing 
the scaling relation pa oc (I-I- 2 )®, constrained to pass through the 
geometric means (dot-dashed lines) of both distributions. 


majority of structures experience an inner density change of 
more than a factor of 2 from 2 = 2 —>■ 0. 

So in this simulated universe, at least, most outcomes 
are possible. Some consistency with our general rules for 
structure formation can be seen, but they are indeed general] 
we should not be surprised when any one particular bit of 
universe - real or simulated - chooses not play by them. 

This brings us back to the particular bits of the simula¬ 
tion that constitute our subsample of compacts. As might be 
expected from the original selection criteria, they are all in 
the upper half of the main panel, all experiencing an increase 
in inner density between the two epochs. This increase is in 
most cases very significant, averaging over a factor of 10 for 
the subsample as a whole, and in several cases exceeding a 
factor of 100. 

Fig. 6 also confirms the feature identified qualitatively 
in Fig. 4, that most of the subsample change relatively lit¬ 
tle in mass ((log[Mv (2 = 0)/Mv{z = 2)]) = 0.12) and some 
30% actually decrease. At the same time, about 10% of the 
subsample exceed the average mass increase for the popula¬ 
tion, thus reinforcing again the point that we should expect 
exceptions. 


■* The mean density change is a little above zero, due to the skew 
in the distribution from systems becoming substructures, but the 
mode is still close to zero. 
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Figure 6. The correlation between mass growth and inner den¬ 
sity change, as seen in structures from the BOLSHOI simulation. 
Background shading in the main panel shows the masses and in¬ 
ner densities of all substructures at 2 : = 0 as fractions of the 
same quantities for their main progenitors at 2 = 2. The dot- 
dashed lines show the overall mean of each quantity. The solid 
and dashed lines show the mean and standard deviation of the 
density change in each mass change bin. Our category of compact 
structures are highlighted as circles, and as dots if they end up as 
substructures at 2 = 2. The peripheral panels show the projection 
of the distributions on to the respective axes, with dotted lines 
showing the distribution for those which become substructures 
by 2 = 0. The white dashed lines and error bars show the mean 
and standard deviation for the compact subsample. The shaded 
area in the right panel highlights the fact that only a minority of 
structures retain their 2 = 2 inner density to within a factor of 2. 

4 SUMMARY 

This study begins by taking the parameters of a known mas¬ 
sive, compact galaxy and looking for simulated structures 
in the BOLSHOI CDM cosmological simulation that would 
broadly correspond with these properties. We Hnd that the 
best-fitting inner density and NFW scale radius of the real 
system would indeed put it amongst the most compact 1% 
of structures in the simulated population from the relevant 
mass-range (Fig. 1). 

Of these analogue compact structures in the simulation, 
the great majority ( A 80%) are substructures within larger 
collapsed regions at 2 = 0, and the majority of the rest were 
substructures in the past. The hosts at 2 = 0 are evenly dis¬ 
tributed in mass, meaning - due to decreasing numbers over¬ 
all - that the probability of finding such a compact structure 
in a given host increases with the host’s mass (Fig. 2) 

In terms of the probability of a substructure of this mass 
being compact - which is perhaps a more relevant statistic - 
we find that this also increases with host mass, the mean ris¬ 
ing from 5 —>■ 30% as Mhost increases from 10^^ —f 10^® M©. 

We then trace the main progenitors of these compact 
analogues back to 2 = 2, where they are all found to be 


distinct systems. At this earlier time, they also had unre¬ 
markable profile parameters (Fig. 4), none of which would 
fall into the extreme category that would later define them 
at 2 = 0. What the subsample did appear to have in common 
was that they are mostly: 

• distinct systems that become substructures, 

• increasing relatively little, or decreasing, in total mass. 

The latter conjecture was quantified, finding that the 
average increase in total virial mass to be just half the aver¬ 
age growth for the population (Fig. 6). Their inner densities, 
meanwhile, had all increased; on average by a factor of 10, 
and many by over 100. 

Concerning the population of simulated structures as 
a whole, we find a clear anti-correlation in the simulation 
between AMv and Aps; mass loss usually leading to higher 
inner densities (more centrally concentrated mass), and late 
mass accretion leading on average to lower eventual inner 
densities (less centrally concentrated). 

In conclusion, the evolutionary paths of simulated struc¬ 
tures revealed even just by this simple analysis, are much 
more rich and varied than our stock phrases on structure 
formation would encapsulate. Even the relatively small vol¬ 
umes of our simulated universes are now big enough that 
improbable things exist, and in great numbers. So we should 
not be afraid to seek them in the real one. 
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mean and standard deviation of this ratio as a function of struc¬ 
ture mass, and also the mean value for substructures. The com¬ 
pact subsample are highlighted as dots. The right panel shows 
the projected distribution in the ratio for the entire mass range. 
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APPENDIX A: CHECKING THE DENSITY 
PROFILES 


The inner densities referred to in all the figures and analysis 
in §§3.1-3.2 are taken from the best-fitting NFW profile to 
each structure: 


p(r) 



(Al) 


The profiles are constrained to satisfy MNFw(fv) = and 
then the best-fitting value of (given this constraint) is 
obtained. The inner density is therefore found from: 

_ Mv __1__ 

47rr| In (1 -|- c) — c/ (1 -|- c) 



For the analysis to carry weight, it is therefore important to 
check whether or not these fits were successful. 

One way of assessing this is to find the radius which 
encloses 2500 times the critical density, r 25 oo, and compute 
the mass predicted to lie within this radius according to the 
best fitting NFW profile. Dividing this expectation by the 
actual mass enclosed, M 2500 , then gives a useful measure of 
the fit accuracy (i.e. a value of one suggests a good fit). 

The result of this simple assessment for the population 
of structures at z = 0 is presented in Fig. Al, showing that 
the enclosed mass predicted by the NFW fit is usually within 
10% of the true value, and almost always within ~ 20%. 
More quantitatively, 98.6% of the structures have a suffi¬ 
ciently good fit to lie within the limits of the y-axis in Fig. 
Al; as can be appreciated by looking at the projected dis¬ 
tribution in the right panel. 

The compact structures we are particularly interested 
in are also representative of this general distribution. The 
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Figure A2. The concentration parameters of simulated struc¬ 
tures as a function of mass, highlighting in particular (shading) 
those whose profiles are poorly matched by an NFW profile, i.e. 
|A^nfw('’250o)/A^ 2500 — 1| > 0.3. As in the other figures, solid 
and dashed lines show the mean and standard deviation of the 
entire sample, and the overall distributions are projected on to 
the right panel. Our subsample of compact systems are also high¬ 
lighted as dots, showing that they are safely at the opposite ex¬ 
treme of the distribution from the poor profile fits. 



Mass growth, M^{z—Q) ! M^(z =2) No. of Structures 


majority of fits to these profiles lead to very nearly the cor¬ 
rect value for M 2500 , and even the few poorer fits are still 
within 20% of the mark. 

Regarding the one or two in every hundred structures 
that have significantly different profiles (i.e. lie outside the 
plot range), these are evenly distributed with mass and there 
is no tendency for being either distinct or substructures. 
What these ‘non-NFW’ systems do have in common is that 
they nearly all produce fits with larger scalelengths, > r^. 

This is shown in Fig. A2, in which it is also interesting 
to see the variation in concentration for the entire sample. 
The projection in the right panel confirms that almost all 
the poorly fitted profiles are low-concentration systems, ap¬ 
proaching or actually hitting® the limit Tg ^ Tv imposed by 
the fitting function. 

The existence of this limit also means that the inner 
densities estimated using these profiles will have a forced 
minimum, slightly greater than the overall virial overdensity: 


P0,min 


1 My 

In 2 — 1/2 drrr® 


1.73 pv(z) , 


(A3) 


This can be seen as the left-hand cutoff in the density dis¬ 
tributions of Figs. 1 and 4 (right-hand panels). The cutoff 
noticeably occurs at lower density at the lower redshift, due 
to evolving Pv(«). 


APPENDIX B: SCALE DENSITY AND 
ENCLOSED DENSITY 

Whether any given aspect of structure growth is best un¬ 
derstood in terms of comoving distances, fixed physical dis¬ 
tances, or some other distance, such as Ts, that reflects the 
self-similarity of the structures, is not always obvious. The 


® About 1000 structures in this My > 10^sample hit this 
limit (i.e. c < 1.1). Of the whole sample, this is barely l/ioooth, 
so no great concern. But of the poor profile fits in general, it 
constitutes about i/ioth. 


Figure Bl. An analogous figure to Fig. 6, using an identical key 
but showing the change within a fixed physical radius for each 
structure, rs [2 = 2] (rather than the change in scale density, as 
before). 


issue can be introduced by rewriting the expression (A2) for 
the inner density of the NFW profile, in terms of Pv, the 
mean density enclosed within r^-. 


_3d^_ 

ln(l c) — c/ {1 + c) 


= p^fic) (Bl) 


i.e. the inner density is a function only of the cosmology and 
of c. So Fig. 6 can also be understood as showing concentra¬ 
tion change as a function of mass growth. 

To help further interpret Fig. 6, we can investigate what 
is happening within a fixed physical radius in each structure 
as it forms, the density enclosed by fixed radius, r, being 
given by: 


p(< r) 


3My /nfw { r / re ) 
47rr3 /nfw(c) 


(B2) 


where /NFw(a:) is the function: ln(l -\- x) — 1/{1 -\- x). The 
change in density enclosed by the original scale length of the 
halo, rs[ 2 ], is therefore: 


Po(< ^^[O]) _ My(0) /NFw(c[a:]) /NFw(ra[2]/rB[0]) 
p7(< My(2) /nfw(c[0]) /nfw(1) 

The distribution of this ratio is shown in Fig. Bl. The key 
result is that mass growth does - on average - not affect 
the inner profiles of the population, over 90% of structures 
preserving the enclosed density within their original scale 
radius to within a factor of 2. 

Our selection criteria, picking out extremely concen¬ 
trated objects at 2 = 0, corresponds here to a subsample 
which do have a slight positive mean fixed density increase. 
But the relative lack of growth in mass - as discussed in §3 
- sets them apart more significantly from the population as 
a whole. 
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